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ABSTRACT 

After a brief review of the theoretical basis of void scaling function properties of 
hierarchical structure, we analyze the phcnomcnological consequences at single jet 
level in Monte Carlo e + e~ annihilation events. We find an interesting alternative 
approach for characterizing quark and gluon jets. 



1. Introduction 

In a field like multiparticle dynamics, where the underlying theory (QCD) 
is as yet too difficult to be useful for predictive calculations at final particle level, 
the study of experimental regularities and their interpretation is of uttermost im- 
portance. As is widely known, regularities which have been discovered at complex 
levels include the NB regularity within final charged hadron multiplicity distribu- 
tions (MD's). More recently, it was found experimentally 1 that separating the 
full sample of events into subsamples with a fixed number of jets, NB regularity 
is satisfied with better accuracy. Our studies of single jets 2 ' 3 shows that the same 
regularity is reproduced at this very elementary level. 

In order to deepen our understanding of single jets MD's, we carried out 
the analysis of the void structure in e + e~ annihilation Monte Carlo events. The 
interesting result is that the differences between quark and gluon jets seen in the 
clan analysis of MD's 2 ' 3 appear relevant also in the 'void' analysis, which becomes 
also an alternative approach to characterizing jets of different origin. 

In Sec. 2 the present theoretical framework for the 'void' function is very 
briefly reviewed; in Sec. 3 our results on Monte Carlo events are commented. 
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2. Hierarchical Correlation Functions and Void Scaling 



The void probability Po(Ay) is defined as the probability of detecting no 
particles in the region of phase space identified by the symbol Ay; here we make 
explicit reference to a central rapidity interval, but it should be clear that any other 
cut can be used: e.g., in Sec. 3 we will discuss the results in rapidity and transverse 
momentum intervals; as another example, in astrophysics Pq is the probability that 
a region of real space is empty of galaxies. 4 

Knowing the void probability is formally equivalent to a knowledge of the 
full MD, since it can easily be shown that 

(-n) n d n 

In practice however it is more useful to consider the void function 4 

V(Ay) = -^—\ogP (Ay) 



ni 

n=l 



where n n (Ay) is the reduced n-order cumulant in the interval Ay defined in the 
standard way as the n-fold integral of the n-order reduced correlation function: 

K n (Ay)= dy 1 ... dy n c n (y x . . .y n ). (3) 

J Ay J Ay 

Notice that there is no average over different intervals of the same size, contrary to 
what is sometimes done experimentally. It will be seen in fact that the regularity 
is better satisfied in central intervals, and degrades when moving away. 

Eq. 3 is very helpful in connection with hierarchical models: these assume 
that reduced correlation functions of any order can be expressed as products of 
two-particle reduced correlation functions, summed over all combinations of pairs, 

c n (y u .. .,y n ) = ^2A n!at ^2c 2 (y il ,yi 2 ) . . .c 2 (yi n _ 1 ,y in ) . (4) 



Here the two-particle reduced correlation functions c^^yi.yj), linking particles of 
rapidity yi and yj, are summed over all non-symmetric relabelings a of the particles, 
and then summed over all distinct topologies a t , with weights A n ^ at which depend 
only on the topology and not on energy or the phase space interval considered. In 
Table 1 the graphical representations of two different hierarchical models are shown: 
edges correspond to two-particle reduced correlation functions, each line corresponds 
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Table 1. The table shows the allowed shapes for a graphical representation of correlation functions 
up to fourth order in the Linked-Pair Ansatz (LPA), and in the Van Hove Ansatz (VHA). 



to a different topology, and different graphs in the same line correspond to different 
relabelings. Two models are shown in Table 1 which differ in the allowed topologies: 
the Linked Pair Ansatz 5 (LPA) requires that a particle appears at most twice in a 
given term of Eq. 4 ('snake' graphs in Table 1), which then becomes simply a sum 
over permutations: 

c n (yi, ■ ■ ■ , 2/tOIlpa =An^2 c 2 (y il ,y i2 ) . . . c 2 (y in _ 1 ,y in ) (5) 

v 

On the other hand, the Van Hove Ansatz 6 (VHA), which was motivated by the 
occurrence of NB regularity, allows in addition connections with three particles 
('star' diagrams in Table 1); in this case one can establish a recurrence relation 
among correlation functions: 

Cn (2/1, ••• ,y n ) I vha = (n - 1) {c n _i (y 1 , . . . , y n -\)c 2 (yt, 2/n)} symm . (6) 

where the right-hand term is symmetrized over all particles. Both these models 
however, coincide at the reduced cumulant level, because integrating Eq. 4 over a 
central rapidity interval Ay one finds the hierarchical structure for cumulants: 

K n {Ay) = A n Kr\Ay) (7) 

It also worth noting that other, non-hierarchical models can approximate relation 
7, thus making it impossible to go backwards from the cumulant level to the corre- 
lation functions level. One of these models, developed in the framework of the 1/H 
expansion, is examined in detail in reference 7. 
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Figure 1. Void function V(Ay) for single jet in e + e~ annihilation (JETSET 7.2) from a) 2-jet 
events, and b) 3-jet events, as a function of n(Ay)K2(Ay) for different central rapidity intervals 
\y\ < 2/cut, with y cu t = 0.25, 0.5, 1.0, 1.5, 2.0 at different cm. energies. The sample of single jets has 
been divided in 5 subsamples as indicated in the figure. Vertical scale refers to the "all sample, 
i.e., the subsample in which no cuts in p t have been applied. Plots for subsequent subsamples are 
shifted down by 0.2. Dashed curves show NB prediction. 



3. Monte Carlo Results 



The considerations developed in the previous section will now be applied 
to e + e~ annihilation events in different rapidity and transverse momentum (p^ 
intervals, with focus on the analysis of single jet properties. 

e + e~ annihilation events have been generated by using Jetset 7.2 (parton 
shower) 8 with Lund string fragmentation as hadronization prescription; the values 
of the parameters of Jetset different from the default ones are listed in reference 2. 
2-, 3- and 4-jet events have been selected by using the Luclus algorithm, and single 
jets have been identified within each sample by the same algorithm; for each sample 
40000 events have been considered at cm. energies y/s = 91 GeV, y/s = 200 GeV, 
y/s = 500 GeV and y/s = 1000 GeV. The analysis has been performed in central 
rapidity intervals \y\ < y cu t with y cut from 0.25 up to the kinematically available 
value and in p intervals p < p with p ^ starting from 0.125 GeV/c. Rapidity 

T T ^Tcut ^Tcut ° I tr J 

and p are defined with respect to the single jet axis. 

In Fig. la we look at the sample of single jets obtained from 2-jet events: each 
event therefore contributes with 2 jets. The figure shows the void scaling function 
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Figure 2. Same as in Fig. 1, but for a) h-jets and b) Z-jets from 3-jet events. In the figure are 
indicated the jet energies Ej corresponding to the cm. energies \fs — 91 GeV (squares), y/s = 200 
GeV (circles), y/1 = 500 GeV (triangles) and yjl = 1000 GeV (diamonds). 

V vs the product nK2- Each point represents the value of the void function for 
particles in a jet of given cm. energy, given rapidity interval and given p T interval. 
For clarity, points referring to different p T intervals have been plotted displaced from 
one another by a fixed amount. The dotted lines, which represent NB behavior, 
become therefore superimposed when this shift is removed. Different center-of-mass 
energies are represented by different symbols, and each point refers to a different 
rapidity interval, up to y cu t = 2.0. We notice firstly that there is a good scaling 
behavior, as the dependence of the function V on energy, rapidity and transverse 
momentum is confined to its dependence on the product n«2- Secondly, we observe 
there is good agreement with the dotted curve which represents NB behavior; this 
was expected on the ground of our previous studies on MD's in restricted domains of 
phase space 2 ' 3 . As for large rapidity intervals, the observed violations to the scaling, 
which are not shown in this figure, are due to the lack of translational invariance for 
the two-particle correlation function which invalidates the hierarchical structure of 
cumulants. Since 2-jet events are composed of a quark and an antiquark jet of the 
same energy, one can conclude that quark jets fulfill NB regularity, and V-scaling 
behavior. 

We now compare Fig. lb, where we show the same quantities as in Fig. la, 
but the single jets here come from 3-jet events. Each event contributes to the sample 



5 



with a quark jet, an antiquark jet and a gluon jet, all of different energies. In this 
figure, all contributions are superimposed. We see that scaling behavior is less good 
than in the 2-jet sample (a fact which also was expected from previous analysis), 
but still acceptable for small p T intervals, and that the data points are spread to 
larger values of fiK2 and at smaller values of V than in the 2-jet sample. Since 
3-jet events contain a gluon jet, we can say that void analysis seems sensitive to the 
presence of a jet originated by a gluon. This is confirmed by the analysis of 4-jet 
events 7 in which the spread is even larger. 

In order to specify better the role of the gluon jet in these events, an energy 
scan of the 3-jet event sample was performed 2 ' 3 . From each event of the 3-jet 
sample the lowest (/-jet), the intermediate and the highest (/i-jet) energy jet have 
been collected in three separated samples. This separation is indeed suggested by 
perturbative QCD: according to conventional wisdom, the behavior of /i-jet samples 
is expected to reflect that of quark (antiquark) jets and the behavior of /-jet samples 
that of gluon jets. The jet energy dependence of the h- and /-jet samples has been 
studied by further collecting the jets in energy intervals 2 GeV wide. The energy 
of each jet is the kinematically reconstructed energy. For this program we analyzed 
100000 3-jet events at cm. energies y/s = 91 GeV, = 200 GeV, = 500 GeV 
and v/i = 1000 GeV. 

Figure 2 shows the results of void analysis on the above sample. In Fig. 2a 
we see the /i-jet sample, and Fig. 2b the /-jet sample. The points are now labeled 
with jet-energy instead of cm. energy, but the same rapidity and p T intervals of the 
other figures have been used. One finds a good scaling behavior and good agreement 
with expectation of NB regularity at all energies, all y cut < 2.0 and all p Tc . It 
should be noticed that these two samples (/i-jet and /-jet) are part of the samples 
that were shown superimposed in Fig. lb. Isolating them and separating them in 
energy has much improved their scaling behavior. In Fig. 3 we compare the two 
samples, the /i-jet and the /-jet at the same jet-energy of 43-45 GeV (of course the 
samples come from events of different cm. energies). The main feature here is the 
different spread of points, because h-jets appear to yield larger values of V than 
/-jets. One should remember that this implies that clans are larger in the case of 
/-jets. 

In order to be assured that /i-jets are actually quark jets, one can compare 
results from the 2-jet sample (where only quark jets appear) and from the h-jet 
sample at the same jet-energy. It is apparent from Fig. 4 that the identification 
of /i-jets with quark jets can be made safely. We have no pure gluon-jet sample 
with which to compare /-jets. However we noticed in Fig. 2b that /-jets also satisfy 
NB regularity, and with different parameters from /i-jets. One expects that if the 
/-jet sample were a mixture of gluon and quark jets, i.e., a mixture of two NB like 
samples with different parameters, NB regularity would not be satisfied. Since on 
the contrary it is, we conclude that in the /-jet sample the contamination of quark 
jets is negligible, and that we can continue to treat /-jets as gluon jets. 
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Figure 3. Comparison of void function V(Ay) 
for h- and Z-jets from 3-jet events as a function 
of n(Aj/)K2(Ay) for the same rapidity intervals 
as in Fig. 1 at jet energy Ej = 43-45 GeV (h- 
jets from cm. energy yfs = 91 GeV, /-jets from 
^fs = 500 GeV). The sample of single jets has 
been divided in 5 subsamples as indicated in 
the figure. Vertical scale refers to the "all p^ 
sample. Plots for subsequent subsamples are 
shifted down by 0.2. Dashed curves show NB 
prediction. 



Figure 4. Comparison of void function V(Ay) 
for single jets from 2-jet events and /i-jets from 
3-jet events as a function of n(Ay)n2(Ay) for 
the same rapidity intervals as in Fig. 1 at jet 
energy Ej — 43-45 GeV (from cm. energy 
\fs = 91 GeV). The sample of single jets has 
been divided in 5 subsamples as indicated in 
the figure. Vertical scale refers to the "all p^ 
sample. Plots for subsequent subsamples are 
shifted down by 0.2. Dashed curves show NB 
prediction. 



We can therefore characterize /i-jets by a little amount of branching which 
grows with energy (clans are small): this situation is consistent with the identi- 
fication of the h-jet as a quark (antiquark) jet within which gluon bremsstrahlung 
is the dominant mechanism. For the /-jet sample branching plays a more relevant 
role than for the /i-jet one: data points spread along the NB curve showing more 
deviation from Poissonian behavior. This result is consistent with the identifica- 
tion of the l-jet as a gluon jet within which gluon self-interaction is the dominant 
mechanism. 

Further support to this result comes from the study of single particle inclu- 
sive rapidity and transverse momentum distributions separately for /i-jets and for 
/-jets; for this program 50000 events have been generated. The inclusive rapidity 
distributions dn/dy are shown vs rapidity y in Fig. 5 at different jet energies: ra- 
pidity is with respect to the ancestor of each jet, event by event; distributions are 
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Figure 5. a) Single particle inclusive rapidity distributions dn/dy vs. y for /i-jets from 3-jet events 
at Ej = 43-45 GcV from y/s = 91 GeV (solid line) and Ej = 248-250 GeV from y/s = 500 GeV 
(dashed line); rapidity is defined in the ancestor frame. Plots are normalized to average multiplicity 
n = 8.8 at cm. energy \fs = 91 GeV and n = 14.5 at y/s — 500 GeV. b) Same as in (a) but for 
/-jets from 3-jet events at Ej = 8-10 GeV from cm. energy yfs = 91 GcV (solid line) and Ej = 
43-45 GeV from = 500 GeV (dashed line). Plots are normalized to average multiplicity n = 7.7 
at = 91 GcV and n = 11.8 at ^ = 500 GeV. 



a b 

Figure 6. a) Single particle normalized inclusive rapidity and p t distributions l/n(d 2 n/ 'dydp^ vs. 
(UiPj) f° r h-jets from 3-jet events at jet energy Ej = 43-45 GeV from cm. energy yfs = 91 GeV; 
rapidity is defined in the cm. frame, b) Same as in (a) but for Z-jets from 3-jet events at jet energy 
Ej = 43-45 GeV from ^ = 500 GeV. 
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Figure 7. Average distance in rapidity between particles for h- and /-jets from 3-jet events as a 
function of jet energy; dotted less dense area shows /i-jets, dotted more dense area shows Z-jets. 

normalized to the average multiplicity. The distribution for /-jets (Fig. 5b) is peaked 
at y = 0, while the distribution for h-jets (Fig. 5a) is more flat. Both distributions 
grow slowly with jet energy. The inclusive (y,p T ) distributions are shown in Fig.s 6a 
(/i-jets) and 6b (/-jets): here the distributions can be compared directly, because 
both have the same jet energy (Ej = 43-45 GeV), both are normalized to 1 and 
rapidity is taken with respect to the cm. frame of the annihilation. The /-jet is seen 
to be less extended in rapidity but more spread out in p , as generally expected for 
a gluon jet with respect to a quark jet. 

In conclusion, all results on single particle inclusive distributions support 
the idea that particles inside a gluon jet, where branching plays a relevant role, 
show more correlation than in a quark jet, where gluon-bremsstrahlung emission 
is dominant. This picture is strengthened finally by Fig. 7, where the average 
distance in rapidity between particles is plotted against jet energy, separately for 
h-jets, which lie all within the less dense area, and for /-jets, which lie all within 
the more dense area, /-jets show more concentration in rapidity (and therefore 
more correlation) than /Vjets, as particles inside a gluon jet are closer together than 
particles inside a quark jet, and in agreement with the fact that the aggregation 
parameter (NB parameter /c -1 ) 9 is larger for /-jets than for /Vjets. 
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4. Conclusions 



The analysis of the structure of voids has been presented as a tool to explore 
the hierarchical structure of correlations in multiparticle dynamics. Its application 
to samples of single jets obtained from e + e~ annihilation Monte Carlo events has 
shown that quark and gluon jets can be characterized by means of the respective 
void properties. It should be noticed that the results obtained are consistent with 
the presence of NB regularity at single jet level, and suggest a hierarchical structure 
for correlations within single jets. 
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